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ABSTRACT 

A recent reanalysis of archival data has lead several authors to arrive at strikingly different 
conclusions for a number of planet-hosting candidate stars. In particular, some radial veloci¬ 
ties measured using FEROS spectra have been shown to be inaccurate, throwing some doubt 
on the validity of a number of planet detections. Motivated by these results, we have begun 
the Reanalysis of Archival FEROS specTra (RAFT) program and here we discuss the first 
results from this work. We have reanalyzed FEROS data for the stars HD 11977, HD 47536, 
HD 70573, HD 110014 and HD 122430, all of which are claimed to have at least one plane¬ 
tary companion. We have reduced the raw data and computed the radial velocity variations of 
these stars, achieving a long-term precision of ~ 10ms' on the known stable star t Ceti, and 
in good agreement with the residuals to our fits. We confirm the existence of planets around 
HD 11977, HD 47536 and HD 110014, but with different orbital parameters than those pre¬ 
viously published. In addition, we found no evidence of the second planet candidate around 
HD 47536, nor any companions orbiting HD 122430 and HD 70573. Finally, we report the 
discovery of a second planet around HD 110014, with a minimum mass of 3.1 Mj^p and a 
orbital period of 130 days. Analysis of activity indicators allow us to confirm the reality of 
our results and also to measure the impact of magnetic activity on our radial velocity measure¬ 
ments. These results confirm that very metal-poor stars down to [Fe/H]—0.7 dex, can indeed 
form giant planets given the right conditions. 
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1 INTRODUCTION 


During the past two decades, the field of exoplanet studies has he- 
come one of the dominant areas of astronomy, leading to the dis¬ 
covery of more 1800 planetfl Out of the methods used to detect 
these systems, one of the most successful for detecting planets or¬ 
biting the nearest stars is the radial velocity (RV) method, which is 
performed by measuring the Doppler shift of the spectral lines of a 
star, produced by the gravitational interaction with orbiting planets. 
The improvement of instrumentation, calibration techniques, and 
signal detection methods has made it possible to detect RV signals 
that conform to planets with masses approaching that of the Earth 
orbiting the nearest stars I Pumusq ue et alfcOlTlnJenkins & Tuomil 
l2014l:lAnglada-Escude et al .112014 . 

iMuller et~ ( 2013 ) has shown that the pipeline of the FEROS 
spectrograph performs an inaccurate barycentric velocity correc¬ 
tion, since it does not include the Earth’s precession, introducing 
a one-year period signal with an amplitude of ~ 62 m s“' for t 


* http://exoplanet.eu 


Ceti (a known stable star at the few m s“' level, see iTuomi et al.l 
l2013h . In addition, the pipeline uses the initial time of the observa¬ 
tion instead of the central time, which also introduces another un¬ 
certainty, especial ly for long exposure ob servations. This issue was 
also mentioned bv ISetiawan et aLloOPOh . This has led some groups 
to question the detections that were made using FEROS data before 
this problem was discovered, leading to the reanalysis of data taken 
with this instrument and the rej ection of some sys tems (e.g. HIP 
13044 in 
|2013| and 


Muller et al.l20ll) . 


_ rei ' 

.Tones & .Tenki^l2014 and HIP 11952 in lPesidera et al.l 


In this paper we present a reanalysis of FEROS data for the 
stars HD 11977, HD 47536, HD 70573, HD 110014 and HD 
122430, all of which were observed with FEROS and are claimed to 
have at least one planet orbiting them. Four of these stars are G and 
K giants and one of them is a very young and nearby star, making 
them all interesting cases in the study of formation and evolution 
of planetary systems. 
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2 OBSERVATIONS AND DATA REDUCTION 

The data were obtained usin g the Fiber-fed Ext ended Range Op¬ 
tical Spectrograph, FEROS Jkaufer et alJll999l) . mounted on the 
2.2m MPG/ESO telescope, at La Silla Observatory. The raw data 
are available in the ESO archive^ The extraction of the FEROS 
spectra was done with the ESO Data Reduction System (DRS), 
available for FEROS users. The DRS performs a bias subtraction, 
flat fielding, order tracing and extraction. The wavelength calibra¬ 
tion was computed using two calibration lamps (one ThAr and 
one ThArNe), having different exposure times, allowing sufficient 
line coverage across all of the spectral range (~3500 - 9200 A). 
The wavelength solution leads to a typical RMS of 0.005 A from 
~900 emission lines. Additionally, the reduction pipeline applies 
a barycentric correction to the extracted data, but this option was 
disabled and applied later using our own code. 

2.1 Radial velocity calculation 

The RVs were c omputed in a similar way as described in 
Ijones et ^ ( l2013h . according to the following procedure. First, 
each order was separated into four chunks an d the doppler shift was 
computed by applying a cross correlation jTonrv & DavisI 0*97911 
between the stellar spectrum and its corresponding template (high 
S/N spectr um of the same sta r). For this purpose we used the IRAF 
fxcor task ( lFitzDatri^ll993h . We only used 35 of the 39 orders 
available, rejecting one in the reddest part of the spectrum (because 
of the presence of telluric lines) and three in the bluest part of 
the spectrum (low S/N). That left us with a total of 140 chunks, 
and therefore 140 different velocities per observation. Then, for 
each RV dataset, the mean velocity was computed, rejecting ev¬ 
ery point more than 2.5 sigma from the mean. The RVs that were 
rejected constitute, on average, 17% of the total number of chunks 
for each epoch. These rejected chunks were the same for most of 
the datasets, at the very blue and red ends of the considered spec¬ 
trum, therefore for most of the epochs we used the same chunks 
when computing the RVs. The principal cause for rejection of a 
chunk was because of the low S/N. 

Similarly, we computed the cross correlation between the si¬ 
multaneous calibration lamp (sky fiber) and one of the lamps that 
was used for the wavelength calibration of that night. The veloc¬ 
ity we measured corresponds to the nightly drift that is produced 
mainly by small pressure and temperature variations inside the 
spectrograph. Finally, we applied the barycentric correction to the 
measured velocities using the central time of the exposure and the 
actual coordinates of the star. 

To test our method we calculated the RV of 127 spectra of t 
Ceti, taken from 2004 to 2013. The resulting velocities are shown 
in Figure [T] The RMS around zero is ~10 m s~ *, which is higher 
than the ~ 4 m s“* obtained bv ijones et alj i2013h for the same star, 
however the difference is due to the limited data used in their work 
that spanned just over 3 years and were observed after 2007. Prior 
to this date there was no standard calibration plan for FEROS, thus 
the wavelength solution was poorer. 

2.2 Data from other instruments 

In this work we also supplemented our data with spectra taken 
with other instruments, namely HARPS, CORALIE and CHIRON. 

^ http://archive.eso.org/eso/eso_ai'chive_main.html 



Figure 1. Radial measurements for r Ceti from 2004 to 2013. The RMS of 
the data is 10 m s“*. 

The High Accura cy Radial velocity Planets Searcher, HARPS 
dMavor et alj 20031) , is located at the 3.6m ESO tele scope a t La Silla 
Observatory. Using this snectrograph. lMtiller et alj ( l2013h obtained 
a RMS of 6.9 m s“* for t Ceti, proving the accuracy of its measure¬ 
ments. The data was obtained from the ESO Reduced Spectral^ For 
each star with HARPS data, the mean value of the RVs was sub¬ 
tracted for each epoch, obtaining velocities distributed around zero 
m s“'. We did this to later combine the HARPS RVs with data from 
other instruments , which were compu ted using a different template. 

CORALIE dOueloz et al.l |20()^ is at the L2m Swiss tele¬ 
scope located at the La Silla Observatory. In ISegransan et al.ld201(]h 
they show that 40% of the measurements made with this instru¬ 
ment reach a radial velocity accuracy of ~ 6 m s“* or better, and 
90% with an accuracy better than ~ 10 m s“'. We obtained the 
CORAL IE RV’s from the Sy stemic Console, which were men¬ 
tioned in ISetiawan et al.l d2008l . hereafter SOS) but were never pub¬ 
lished. 

Einally, we used data from the CHIRON spectrograph 
dTokovinin et alJ l20I3h . located at the L5m telescope at Cerro 
Tololo Inter-American Observatory (CTIO). The extraction of the 
data was done using the automatic reduction pipeline that is offered 
to CHIRO N users. The veloc ities were obtained using the U cell 
technique dButler et alJll996l) . This consists of passing the stellar 
light through a cell that contains iodine vapor, that superimposes 
absorption lines in the region between 5000 A and 6000 A. These 
lines are then used as markers against w hich the doppler sh ift can 
be measured. More detail can be found in llones et al.ld2014h . With 
this instrument following our procedures it is possible to reach a 
velocity precision of ~ 6 m s“*. 

The fits for the RV cur ves were made using the Systemic Con¬ 
sole dMeschiari et al.ll2009h . We changed the objective function of 
the minimizer to be the RMS of the fit, instead of they^. This was 
done because of the difference between the uncertainties obtained 
in the measurements for each instrument. The HARPS data have 
uncertainties on the order of ~ 0.5 m s“*, whereas for the other 
datasets they are generally higher than 5 m s“*. When we tried 
to minimize the y“, the HARPS data would dominate the fitting 

^ http://archive.eso.org/wdb/wdb/adp/phase3_spectral/form? 
phase3_collection=HARPS 
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Figure 2. Top panel. Radial measurements for HD 11977 taken with 
FEROS (red circles), HARPS (blue squares), and CHIRON (green trian¬ 
gles). The back solid line corresponds to the RV curve found in this paper, 
with a RMS of 11.2 m s“ ^. The dashed grey line corresponds to the fit found 
in S05, which gives us a RMS of 33.8 m s“^ Bottom planet Residuals from 
our fit. 


Table 1. Orbital parameters for HD 11977 b 


P (day) 

621.6 ±2.4 

To (JD-2450000) 

2906.6 ± 10 

e 

0.3 ± 0.03 

OJ (deg) 

35.5 ± 5.8 

Mp sin i (Mjup) 

6.5 ± 0.2 

a(AU) 

1.89 ±0.005 


process, even when most of the RVs were computed from FEROS 
spectra. Minimizing the RMS prevents this from happening, mak¬ 
ing the fit more evenly distributed across all the datapoints. The 
uncertainties in the derived orbital parameters for each star were 
obtained using the boot strapping option available in the Systemic 
Console. 
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Figure 3. LS periodogram of our RV data for HD 11977. The two horizontal 
lines represent the 1% (solid line) and 0.1% (dashed line) FAR 


eludes a total of 48 spectra. We also included 217 measurements 
taken with HARPS. Since many of these RV epochs were obtained 
during a single night, we decided to use the RV average in that 
night, leaving us with 13 RV epochs. We also included 8 observa¬ 
tions made with CHIRON. Figure shows the resulting RVs from 
FEROS (red circles), HARPS data (blue squares), and CHIRON 
(green triangles). When fitting the orbital solution found in S05, we 
got a RMS of 33.8 m s“*, larger than the one published i n S05 of 
29.1 m s '. Figure[3]shows the Lomb-Scargle periodogram dScarglel 
Il982l) of the data. As can be seen, there is a strong peak at 625 days. 
After close inspection, it was found that this period does not corre¬ 
spond to any peak in the periodogram of the sampling, however we 
did find that the second peak at 23 days is a product of the sampling 
of the data. We fitted a single-planet Keplerian solution to our data, 
starting from the 625 day period found in the periodogram. The 
best fit leads to the following parameters: P = 621 days, MpSini 
= 6.5 Mjup, and an eccentricity of e = 0.3. The orbital parameters 
are listed in Table[T] This solution leads to a RMS of 11.2 m s ', 
almost a factor three better than the 29.1 m s“' uncertainty found 
in S05. 


3 RADIAL VELOCITY ANALYSIS 
3.1 HD 11977 

In ISetiawan et al.l ( 120051 hereafter S05), this star was classified as 
a G5III star, located in the ’’clump region”. They derived a mass of 
1.91 Mq, a rotational period of Pj-Qt/sini = 230-270 days, and a 
metallicity of [Fe/H] = -0.21. Assuming negligible mass loss dur¬ 
ing the red giant branch phase, it is likely a former main-sequence 
A-star. We will adopt a mass of 2.31 Mp, and metal licity of [Fe/H] = 
-0.16, which were derived in iMortier et alj ^20131 hereafter M13). 
S05 also announced a planetary companion with an orbital pe¬ 
riod of P = 711 days, an eccentricity e = 0.4, and minimum mass 
m 2 sin i =6.5 Mjup, being the first planet found orbiting an interme¬ 
diate mass star (for M* = 1.19Mo). 

To test this Keplerian solution we used FEROS data of this star 
taken between 2001 and 2007. However, when testing the nights 
that were observed before 2004 we obtained a large scatter in the 
velocities for t Ceti (~ 30 m s“*), making those nights unreliable. 
Therefore, we only used the nights from 2004 to 2007, which in- 


3.2 HD 47536 

The se cond star in our sample is HD 47536, a very metal-poor Kl- 
III star. lSetiawan et all l l2003l) derived a mass of 1.1-3.0 Mq, a rota¬ 
tional period of Prot/ sin i = 619 days, and a metallicity of [Pe/H] = 
-0.68 dex. It was first studied using data from 1999 to 2002 taken 
with EEROS, and one night with CORALIE in October 2002. They 
reported the discovery of one companion with a period of 712 days 
and a minimum mass of 5.0 Mj^p. Later, S08 did a revision of the 
orbital solution including new data from 2004. They claimed that 
the star actually hosts two planets with the following orbital pa¬ 
rameters: P ~430 days and Mp ~5 Mj^p for the inner planet, and a 
P ~2500 days and a Mp ~7 Mj^p for the outer object. It is worth 
mentioning that they did not publish the exact results from the two- 
planet solution, the values we quote came from an unrefereed con¬ 
ference proceeding. 

We analyzed 56 EEROS spectra of HD 47536, taken from 
2005 to 2011. We did not include the data taken before 2005 be¬ 
cause of the large RV scatter we obtained for t Ceti, as was the 
case for the nights before 2004 for HD 11977. Additionally, we 
included 12 RV epochs obtained with CORALIE, 6 data points ob- 
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Figure 4. Top panel. Radial velocity measurements for HD 47536 
from FEROS (red circles), CORALIE (orange diamonds), HARPS (blue 
squares), and CHIRON (green triangles). The back solid line corresponds 
to the RV curve found in this paper. The RMS of this fit is 51.7 ms“^. The 
dashed grey line corresponds to the fit found in ISetiawan et alj j2QQ3h . with 
a RMS of 103.5 m s“^. Bottom panel'. Residuals from the fit found in this 
paper. 




Figure 5. Lomb-Scargle periodogram of the HD 47536 RV measurements. 
The two horizontal lines represent the 1% (solid line) and 0.1% (dashed 
line) FAR 


tained with CHIRON, and 18 RVs computed from HARPS spectra. 
The resulting RV measurements are shown in Figure [4] The red 
circles, orange diamonds, blue squares, and green triangles corre¬ 
spond to FEROS, CORALIE, HARPS, and CHIRON data, respec¬ 
tively. We also used the mass and metallicity derived by M13 for 
this star, corresponding to M* = 0.98Mq and [Fe/H]=-0.65 dex. 
The LS periodogram of the RVs (Figure[3 shows a strong peak at 
434 days. Starting from this period we fit a Keplerian solution to the 
data, which leads us to the following orbital parameters: P=434.9 
days, e = 0.3 and M/> = 4.0 Mjup. The orbital solution is overplotted 
in Figure|4](solid line) and we obtain a RMS around the fit of 51.7 
ms“^ The orbital parameters are listed in Table 

In Figure[5|the periodogram for window function and the data 
are shown. We see that there is a peak in the sampling, at 454 days, 
very close to the peak found in Figure [5] In order to discard the 


Table 2. Orbital parameters for HD 47536 b 


P (day) 

434.9 ± 2.6 

To (JD-2450000) 

3040.8 ± 30 

e 

0.3 ±0.1 

(D (deg) 

268.6 ± 17.3 

Mp sin i (Mjup) 

4.0 ± 0.4 

a(AU) 

1.12 ±0.005 


434.78 d 
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- 
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Figure 6. Lomb-Scargle periodogram of the sampling window (black line) 
and the data (grey line) for HD 47536. The peak at 434.78 days corresponds 
to the peak found in Figure!^ 


possibility that the 434-days period is not a product of the sam¬ 
pling of the data, we measured its power in the periodogram as we 
add data points. If the period is from a true signal from the star, 
then we expect that it should become stronger as we add more RV 
measurements, in a monotonic fashion. If it is a product of the sam¬ 
pling, its power could stop growing and will decline in regions. We 
found that the power grows as we add more data points more or less 
monotonically, instead of decreasing, which helps to confirm that 
the 434-days period is produced by a true signal from the system. 

We then tried to fit a two-planet solution using as initial 
guesses the values for the period and mass of the planets found 
in SOS (limiting the eccentricity of the outer planet to he smaller 
than 0.5 otherwise the system might be highly unstable). We found 
a mass for the outer planet of ~ 2.0 Mj^p, with a period of 2013 days 
and eccentricity of 0.5. This two-planet fit is shown in Figure|7] The 
RMS of the solution is 48.9 m s“*, which is slightly better than the 
one-planet fit. We performed an F-test between these two solutions 
and found that the two-planet solution is not statistically superior 
to the one-planet fit (probability of 70 %). 

When looking again at the periodogram of the data, we were 
not able to find any other significant peak that could justify another 
signal. The periodogram of the residuals from the one-planet fit 
(Figure [8) does not show a significant peak neither. Therefore, we 
can rule out the existence of a second statistically significant signal 
in the current RV data of this star. 


3.3 HD 110014 

HD 110014 is a K2III star with an absolute magnitude = -0.11, 
color [B - V) = 1 .24, and metallicity [Fe/H]=0.14 dex (Ml3). In 
Ide Medeiros et al.l (120091 hereafter M09) they estimated the upper 
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Figure 7. Top panel. Radial velocity for HD 47536 from FEROS (red cir¬ 
cles), CORALIE (orange diamonds), HARPS (blue squares), and CHIRON 
(green triangles). The black solid line corresponds to the RV curve found 
in this paper for a 2-planet fit, with a RMS of 48.9 m s“^ The dashed grey 
line corresponds to the fit found in SOS, with a RMS of 58.8 m s“^. Bottom 
planet Residuals from the fit found in this paper. 
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Figure 8. Periodogram of the residuals of HD 47536 after fitting the one- 
planet solution. The two horizontal lines represent the 1% (solid line) and 
0.1% (dashed line) FAR 


limit of the rotational period as Prot/ sini ~ 513 days. They also 
estimated the mass and age of this star by its positio n in the color- 
magni tude diagram, using the process explained in Ida Silva et al.l 
hoOhl) . According to their solution, HD 110014 could be either a 
1.9 Mo star at the end of its core helium phase, or a 2.4 Mo star in its 
first-ascent of the RGB, close to the core helium flash. In this paper 
we will adopt the mass of 2.09 Mo, as derived by M13. In M09 
this star was observed using FEROS from 1999 to 2007, plus a few 
nights with HARPS and CORALIE. They claim the detection of a 
planet with a period of 835.5 days, a minimum mass of 11.1-9.5 
Mjup (for M* = 1.9 - 2.4 Mq), and an eccentricity of 0.46. 

We analyzed spectra taken with FEROS from 2004 to 2011, 
corresponding to 21 nights and a total of 25 data points. We also 
used 116 data points taken with HARPS, but as with HD 11977, 
many of them were taken in the same night and so we took the 
average velocity per night. That left us with 17 HARPS data points. 



Figure 9. LS periodogram of HD 110014. The two horizontal lines repre¬ 
sent the 1% (solid line) and 0.1% (dashed line) FAR. 



Figure 10. Top panel: Radial measurements for HD 110014 from FEROS 
(red circles) and HARPS (blue squares). The back solid line corresponds to 
the RV curve found in this paper, with a RMS of 44.6 m s“*. The dashed 
grey line correspond to the fit found in M09, which gives us a RMS of 63.13 
m s“'. Bottom planet: Residuals from the fit. 


We tried to fit the solution found in M09 to our data, obtaining 
a RMS of 63.13 m s“*, larger than the one published of only 45.8 
m . When trying to find another orbital solution we looked at the 
periodogram of the data in Figurej^ which shows an emerging peak 
at 833 days. Beginning with that initial period and minimizing all 
the orbital parameters, we obtained a fit with a RMS of 44.64 m s“*. 
The fit is shown in FigurefTOland the orbital parameters are listed in 
Table[3 The period of this possible companion is 936 days, slightly 
larger than the 833-days period found in Figure|9] The periodogram 
of the residuals (Figure fTTt shows a possible peak at 133 days and is 
comfortably above the 0.1% FAR When a new signal with a period 
of 133 days is added to our previous one-planet solution, we find 
a new model as shown in Figure [T2l and this leads to the orbital 
parameters listed in Table |4] The RMS of this two-planet solution 
is 19.4 m s“*, significantly lower than the one found in M09 of 45.8 
m s“' and our own one-planet solution. 

We performed a F-test between the one-planet solution and 
the two-planet solution and found that they are significally different 
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Table 3. Orbital parameters for HD 110014 b. 


Table 4. Orbital parameters for HD 110014 b and HD 110014 c. 


P (day) 

To (JD-2450000) 
e 

oj (deg) 

Mp sin i (Mjup) 
a(AU) 


936.4 ± 44.5 
3842.0 ± 206.6 
0.06 ± 0.09 
119.55 ± 84.2 
13.67 ± 2.3 
2.5 ± 0.08 



b 

c 

P (day) 

882.6 ±21.5 

130.0 ±0.9 

Tq (JD-2450000) 

3739.4 ± 50.2 

4031.3 ±3.78 

e 

0.26 ±0.1 

0.44 ± 0.2 

at (deg) 

47.72 ± 15.5 

307.6 ± 19.9 

Mp sin i (Mjup) 

10.7 ± 1.0 

3.1 ±0.4 

a(AU) 

2.31 ±0.04 

0.64 ± 0.003 



1 10 100 1000 10000 
Period (days) 


Figure 11. Lomb-Scargle periodogram of the residuals of HD 110014, after 
fitting the one-planet solution. The two horizontal lines represent the 1% 
(solid line) and 0.1% (dashed line) FAR 


(probability of ~0% that they are both statistically similar). We also 
looked at the periodogram of the sampling (Figure fOl l and we find 
that there is no peak corresponding to the 133-days period. This 
helps us to rule out the possibility that the signal we see is due to 
the sampling of the data. 

In M09 a periodicity of 130-days was also detected in the 
residuals of their fit. They performed several tests in order to dis¬ 
card the possibility that the signal is intrinsic to the star. They could 
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Figure 13. Lomb-Scargle periodogram of the sampling of the data for HD 
110014. The vertical dashed lines represent the location of the 833-days 
period and the 133-days period found in the periodograms of Figures 
and l 11 1 


not find any correlation between those tests and the measured ve¬ 
locities, nor with the fit residuals, yet, they could not find enough 
evidence to he sure that the signal is produced hy the presence of 
a second companion. The RMS they obtain after fitting the orbital 
solution of the second planet is 32.9 m s“*, significantly higher than 
the 19.4 m s“* RMS we got with our fit. 

The 130-days period could be caused, not only by a planetary 
companion, but also by a cool spot on the surface of the star. It has 
been shown that RV signals caused by rotational modulation of star 
spots can mimic Doppler sig nals, causing possible misinterpreta¬ 
tion o f the measurements (e.g. lSetiawan et al.l2008l ; lHuelamo et al.l 
l2008h . particularly when the star is at low inclination since some 
activity indicators, li ke the bisector inv erse slope (BIS), are insen¬ 
sitive in this regime dPesort et al.ll2007l) . In order t o test this possi - 
bility, we first computed the spot filling factor from lHatz^ ||2002|) . 
This is the percentage of the surface of the star that a spot should 
have to cover in order to produce a RV signal with the same am¬ 
plitude as the one measured. We us ed V = vsini = 1.1 km s“*, 
derived in IPe Medeiros et all ( l2014h . as the rotational velocity of 
the star, and A '2 = 83.5 m s“* as the RV amplitude produced by the 
130-days signal. We obtained a filling factor of 13.8%. In section 
14.21 we discuss the analysis of available photometry for this star, to 
test if we can find signals that could be due to spots with this level 
of spot coverage. 


Figure 12. Top paneT. Radial velocity measurements for HD 110014 from 
FEROS (red circles) and HARPS (blue squares). The back solid line corre¬ 
sponds to the RV curve found for a two-planets system. The RMS of this fit 
is 19.4 m s“*. Bottom planet Residuals from the fit. 


3.4 HD 122430 

HD 122430 is a giant star (K2III) with metallicity of [Fe/H] = -0.09 
dex and a mass of 1.68 Mq (data from M13). A planet orbiting this 
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Figure 14. Lomb-Scargle periodogram of HD 122430. The two horizontal 
lines represent the 1 % (solid line) and 0.1% (dashed line) FAR 



Figure 15. Top panel: Radial velocity measurements for HD 122430 from 
FEROS (red circles) and HARPS (blue squares). The back solid line corre¬ 
sponds to the RV curve found in this paper a nd has a RMS of 29.3 m s“'. 
The grey dashed line is the solution found in lSetiaw^ (200-3h and gives a 
RMS of 63.4 m s“*. Bottom planet Residuals from the fit. 



Figure 16. Top panel: Radial velocities measured by FEROS for HD 70573 
(red dots). The grey line correspond to the orbital solution found in S07. 
The RMS of this fit is 78.8 m s“'. Bottom panel: Residuals from the fit. A 
long period variation can be seen in the residuals. 


still lower than the one obtained for HD 47536 and HD 110014 
(one-planet fit), their periodograms showed at least one period with 
a significant peak close to the 1% and 0.1% FAP, which helped us 
prove that there was a signal within the RV data. This is not the case 
for HD 122430, and after finding a RMS of 29 m s“* for our best 
solution, we cannot support the existence of this planet based on 
our data. Indeed, a mean noise model fit to the data returns an RMS 
of 40 m s“*, only ~ 10 m s“* larger than the one-planet solution. 

In order to check if we missed the signal in the 
LS periodogram because of its high eccentricity, we also 
calculated the Generalized Lomb-Scargle Periodogram (GLS, 
IZechmeister & Kursteill20()^ . This periodogram is more suited to 
detecting signals arising from bodies on eccentric orbits (e > 0.5). 
We found no peak above the noise level, meaning that the high ec¬ 
centricity was not the only reason no peak was found in the LS 
periodogram. Therefore, if a planet does exist orbiting this star, we 
do not have the necessary data yet to make this claim. 


3.5 HD 70573 


star was reported in lSetiaw^ ll2003h . using data from FEROS and 
CORALIE. The orbital parameters correspond to a planet with a pe¬ 
riod of 344 days, eccentricity of 0.68, and minimum mass Mp sin i 
- 3.71-6.04 Mjup (for M* = L2-2.5 Mo). The RMS of this solution 
is reported to be 27.2 m s“*. 

We used spectra taken with FEROS from 2004 to 2010, to¬ 
taling 42 data points. We also used 6 velocities computed with 
HARPS as a check on our soluti on. When trying to fit the values for 
the orbital parameters found bv ISetiaw^ lioOsL hereafter S03) we 
found a RMS of 63.4 m s“*. From the periodogram (Figure fT4ll we 
could not find any significant peak that would lead us to a starting 
period for a fit. We then let the console minimize the orbital param¬ 
eters with starting values given by the orbital parameters found in 
S03. This leads us to a new solution with a period of 455.1 days, 
an eccentricity of 0.6, a minimum mass of 2.1 Mj^p and a RMS of 
29.3 m s“*. The fit is shown in Figure [TS] as the black solid line, 
where the red circles and the blue squares represent the FEROS 
and HARPS data, respectively. Even though the RMS of S03 is 


The last star in our sample is HD 70573, a young lithium bear¬ 
ing star classified as a member of the L ocal Association (Pleiades 
moving group) by iMontes et alJ l l200lh . The age is therefore es¬ 
timated to be between 78-125 Myr (only 3%-6% the age of the 
Sun, whe re moving groups ar e thought to be associations of young 
stars, see lMurgas et al.|[2013h . It h as a spectral type of G 1-1.5 V 
and a mass of 1.0 Mp, according to lSetiawan et al.1 (120071 hereafter 
S07). lHenrv. Eekel & Halil ( Il995h measured photometric variations 
of HD 70573 and found a rotational period of 3.296 days. 

This star was studied in S07, in which they claimed the de¬ 
tection of a companion with a period of 852 days and a minimum 
mass (Mp sin i) of 6.1 Mj^p. To verify this detection we analyzed 
data taken with EEROS between 2004 and 2008, corresponding to 
a total of 41 nights and 55 data points. Since this is a young star 
some parts of its spectrum contain very few lines or heavily blended 
lines (because of line broadening from rapid rotation). This leads to 
an innacurate correlation between the spectra and the template for 
some orders, and in turn in noisy results for the velocities. To deal 
with this, we had to modify the orders that we were considering to 
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Figure 17. LS periodogram of HD 70573. The two horizontal lines repre¬ 
sent the 1% (solid line) and 0.1% (dashed line) FAR Only one significant 
peak, above the 0.1% FAR, is seen, corresponding to 3.296 days. 


833.33 d 



Period (days) 


Figure 19. LS periodogram of the sampling window for HD 70573. The 
peak corresponding to 833 days, close to the period of 852-days found by 
S07, is shown. 



Phase 


Figure 18. Top panel: Phase-folded RV curve for HD 70573 (red dots), and 
fit to the data (solid line).The RMS of this fit is 36.4 m s“*. Bottom panel: 
Residuals from the fit. 


calculate the RVs. After visually inspecting each order for different 
observations of this star, we decided to consider only from order 9 
to 31, instead of using information from order 2 to 36 like we did 
for the above stars. 

When we fit a companion with the orbital parameters as given 
by S07 to our data we arrived at a RMS of 78.78 m s“', much higher 
than the published RMS of 18.7 m s“*. The fit is shown in FigurefT^ 
as the grey solid line, and the red circles represent the FEROS RV’s. 
The periodogram of the data (Figure [T tIi does not show any signif¬ 
icant peak at high periods, like it did in S07. Only one peak above 
the 0.1% FAP can be seen, corresponding to a period o f 3.296 days, 
the sa me as the rotational period of the star found in iHenrv et aP 
( Il995h . We tried to fit a curve with a period of 3.296 days to model 
the rotational effects of the star on the RVs. The fit is shown in Fig- 
ure[T^and has an associated RMS of 36.4 m s“*. The periodogram 
of the sampling (Figure [T^ shows that one of the peaks correspond 
to a 833-day period, close to the 852-days period proposed by S07. 
This indicates that the previously published orbital solution might 
be due to a sampling alias combined with noisy data. 


4 ACTIVITY 

4.1 S-index and Bisector Velocity Span 


It has been shown that intrinsic stellar phenomena such as mag¬ 
netic activity can produce a periodic radial vel ocity signal that ca n 
be misinterpretated as a substellar companion l lOueloz et alj200lh . 
We performed three tests of our data to verify if this is the case for 
the stars we are studying. First, we measured the change in the flux 
from the Ca II HK lines over time by calculating the S-index fo r 
each epoch (the process is described in ijenkins et al.ll2008L 2011). 
The s econd test we used was the bisector analysis JToner & Gravl 
Il988l) . aimed at detecting asymmetries in the line profiles caused 
by intrinsic stellar phenomena. We measured the bisector velocity 
span (BVS) for each epoch, which essentially corresponds to the 
velocity difference befween th e boftom and top thirds of the cross 
correlation function (CCF, see ljones et al.ll2013h . Finally, the third 
test we used to parameterize line asymmetries was analysis of the 
CCF full width at half maximum (FWHM). 

Our aim was to test if these activity indicators were correlated 
with the RV at each epoch. To measure the level of the correla¬ 
tion we computed the Pearson rank correlation coefficient, which 
is a measure of the linear correlation between two variables, where 
1 is a total positive correlation, 0 is no correlation, and -1 is total 
negative correlation. When computing the Pearson correlation co¬ 
efficient all the data points have the same weight, meaning that one 
point lying very far away from the mean is going to affect the final 
result, even if it has a large uncertainty. That is why we only used 
the points that were closer than 2.5cr from the mean value of the 
distribution. 

For most of the stars no correlation was found between the ac¬ 
tivity indicators and the velocities, with r ~ 0, where r is the Pear¬ 
son correlation coefficient. Plots of the indices vs RV for each star 
can be found in the Appendix. There were some cases that required 
further analysis. One of them is HD 47536. The RVs obtained for 
this star showed a lot of jitter, even after fitting a one-planet so¬ 
lution (~ 50 m s“'), leaving us with a RMS of 51.7 m s“* which 
could not be significantly improved by adding a new component 
to the system. In the left panels of Figure we show the activity 
indicators versus the measured RV, and in the right panels the same 
activity indicators but this time versus the residuals left over after 
fitting our one-planet solution. There are no significant correlations 
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Figure 20. S-index (upper panel), bisector velocity span (middle panel), and 
FWHM of the CCF (lower panel) versus the measured RV (left panel) and 
versus the residual RV for HD 47536, after fitting the one-planet solution. 
Only the values for the activity indices laying 2.5cr from the mean value of 
the distribution were considered. The correlation coefficients for the plots 
in the left panels are r= -0.3, 0.3 and -0.02, and for the right panels are r = 
-0.2, 0.24 and -0.23 for the S-index, BVS and FWHM versus the residual 
velocity, respectively. Horizontal dashed lines correspond to the mean of 
the distribution for each plot. 


Figure 21. S-index (upper panel), bisector velocity span (middle panel), 
and FWHM of the CCF (lower panel) versus the measured RV (left panel), 
and versus residual RV (right pannel) for HD 110014, after fitting the two- 
planet solution. Only the values for the activity indices laying 2.5cr from the 
mean value of the distribution were considered. The correlation coefficients 
for the plots in the left panel are r=0.36, 0.19 and -0.17, and in the right 
panel are r = 0.3, -0.18 and 0.68 for the S-index, BVS and FWHM versus 
the velocity, respectively. Horizontal dashed lines correspond to the mean 
of the distribution for each plot. 


in any of the plots versus the RVs, with r indices equal to -0.3, 0.3 
and -0.02. From the;^^ test between the indices and a mean noise 
model fit, we found that the probability of these being fit by a hori¬ 
zontal line is ~ 90%. This means that it is unlikely that the veloci¬ 
ties are a product of the stellar activity, and are likely produced by 
an external source. Also, we find no correlation in any of the plots 
versus the residual RV (right panel), which might indicate that the 
jitter in the residuals is not due to magnetic activity nor rotationally 
modulated spots on the star. We computed th e velocity oscillation 
amplit ude for this star from the relations in iKieldsen & Bedding 
(Il995h. We used M = 0.98 Mq from M13 and L = 175.19 Lq from 
I Anderson & Fran^(l2012h . We obtained Vosc = (41.8 ±2.5) m s“*. 
This is only 10 m s“' lower than our RMS of 51.7 m s“', meaning 
that the residuals from our 1-planet fit could be explained by the 
expected "jitter” for this giant star. 

In Figure[2T|we can see the same plots as before, but this time 
for the star HD 110014. In the left panels we show the activity 
indices versus the measured RV, and in the right panels are the same 
activity indices but versus the residual velocity left after fitting the 
two-planet solution we are proposing in this paper. We found that 
the r values for the activity indices and velocity in the left panels 
plots are r=0.36, 0.19 and -0.17 for the S-index, BVS and FWHM, 
respectively. This indicates that activity in the star is not the cause 
of the velocities we observe and supports our argument that they 
are produced by two planets orbiting the star. In the right panels, 
the r indices are r = 0.3, -0.18 and 0.68 for the S-index, BVS 
and FWHM versus the residual velocity, respectively, suggesting 
that part of the jitter we were left with after fitting our two-planet 
solution may be caused by activity in the atmosphere of the star. 

The last case is for HD 70573. One interesting result is that, 
when we studied the periodogram of the BVS (Figure 122b . we 
found a signal corresponding to the orbital period of the star of 
3.296 days, just like we did with our RVs. Even though this peak 
is not nearly as significant as it is in Figure [T3 it is present and it 



Period (days) 


Figure 22. LS periodogram of the bisector velocity span for HD 70573. 




Figure 23. BVS vs measured RV (left panel) and vs residual velocity after 
removing the rotational solution (right panel) for the star HD70573. The r 
coefficients are 0.44 (left panel) and 0.26 (right panel). 
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Figure 24. V-magnitude variation for HD 11977 (upper panel), HD 47536 
(middle panel) and HD 110014 (lower panel). In each plot a long-period 
variation was found (red lines), with P = 2114 days for HD 11977, P = 
1962 days for HD 47536, and P = 1687 days for HD 110014. 

supports our result that the observed velocities contain a signal that 
is due to features on the star modulated by its rotation. The S-index 
and FWHM of the CCF do not show correlation with the RV, with r 
indices less than 0.2. Flowever, for the case of the BVS (Figure[23]l, 
we have a correlation coefficient of r = 0.44 for this index vs the 
measured RVs, which is suggestive of correlation. The r value de¬ 
creases to r = 0.26 when considering just the residual velocities 
after fitting the rotational period of the star. This strongly suggests 
that the rotation of the star is likely the cause of the variation seen 
in the bisector. 

4.2 Photometry 

We looked for any period variation in the visual magnitude that 
could mimic the results we show in this work. We anal yzed the V- 
band photometry from the All Sky Automated Survey dPoimanskil 
Il997h for the stars HD 11977, HD 47536 and HD 110014. The pre¬ 
cision from this data is on the order of 1%. The magnitude versus 
time are plotted in Figure [24] Superimposed on each plot are the 
long-period fits, with periods of P = 2114, 1962 and 1687 days 
for HD 11977, HD 47536 and HD 110014, respectively. Because 
these are evolved stars, this long-period variation could be caused 
by magnetic cycles in the photosphere of the stars. These periods 
are not present in the velocities so we can conclude that this photo¬ 
metric variation is not the cause of the signals in the RV data. 

For HD 110014 we analyzed the photometry in order to search 
for any signals from cool spots on surface of the star that could 
cause the RV signal we found (see section [T3l l. We could not find 
any period above the noise level that could correspond to the peri¬ 
ods of th e planets detected. We also looked at the Hipparcos pho¬ 
tometry jPerrvman et al.llT997h for this star but we found no sig¬ 
nificant peak with a period longer than 2 days. The precision of 
the data for this star is ~0.1%. As previously mentioned, the filling 
factor of a spot on the surface would need to be 13.8% to produce 
this level of RV signal, and therefore it is likely that this level of 
spot coverage would not be detectable in this dataset, meaning a 
spot origin for the signal can not be ruled out with any degree of 
certainty. However, we do note that we studied the power of the 
periodogram peak as a function of time, as we did for HD 47536 in 


section|T2] We found that, in general, the signal peak rising almost 
monotonically. This lends some weight to the planetary hypothesis 
for the second signal in the RV data of this star. 


5 CONCLUSIONS AND DISCUSSION 

We computed precise radial velocities using archival FEROS spec¬ 
tra for stars HD 11977, HD 47536, HD 70573, HD 110014 and HD 
122430. From this data we can conclude the following: 

• There is a companion around HD 11977 with a period of 621.6 
days and minimum mass of 6.5 Mj^p. 

• We can confirm the presence of one companion around HD 
47536, with a period of 434.9 days and minimum mass of 4.0 Mj^p. 
We cannot confirm the existence the second companion proposed 
by SOS. 

• We confirm the existence of a companion around HD 110014, 
with a period of 882.6 days and minimum mass of 10.7 Mj^p. We 
also found evidence of a second companion around this star with 
a period of 130.0 days and minimum mass of 3.1 Mj^p. A more 
extensive photometric analysis has to be done to confirm that this 
signal corresponds to a second companion and not to a cold spot on 
the star’s surface. 

• There is no signal significant enough to confirm the presence 
of any companion around HD 122430. 

• We found no evidence of the planet orbiting HD 70573 pro¬ 
posed in S07, nor evidence of any other companion. The only signal 
present corresponds to the rotational period of the star. 

• There is no correlation between the activity indices and the 
RVs for the stars HD 11977, HD 47536, HD 110014 and HD 
122430, supporting our hypothesis that the observed velocity sig¬ 
nals are Doppler signals. We find a small yet not significant corre¬ 
lation between the activity indices and the measured RVs for HD 
70573, but it decreases after removing the rotation of the star. We 
find a higher level of correlation in the activity indices versus the 
residual velocities after removing the two-planet fit for HD 110014. 
This could explain the jitter in the residual RVs as a product of ac¬ 
tivity within the star. 

• We observe a periodic variation in the photometric data for 
three of the stars in our sample (HD 11977, HD 46536 and HD 
110014), but with periods different than the ones of the compan¬ 
ions we found. We conclude that the signals in HD 11977 and HD 
47536 data cannot be explained by intrinsic photometric variation 
in the stars. For HD 110014 we could not find evidence of any pho¬ 
tometric signal with periods equal to the one of the RV signals, 
supporting our argument of the existence of this second planet or¬ 
biting this star. 

Finally, in Figure we can see that the planets found in 
this work follow similar trends to the previously detected planets 
around giant stars. The upper and middle panels show that these 
planets are mostly clustered around periods from 100 to 1000 
days, and with masses from ~ 1 to 30 Mjun. The bottom panel of 
Figure is similar to Figure 2 in iReffert et al.l ( 120151 hereafter 
R15) which shows the distribution of metallicity and stellar mass 
for planets around giant stars. HD 110014 b and c fall in the region 
corresponding to a 20.8% probability of the star having a planet, 
and HD 11977 b falls in the block with probability of 6.9%. Both 
stars are in the range of masses for giant stars for which we can 
expect the highest number of planets (~2 Mq according Figure 2 
and 5 of R15). HD 110014 is also the star with the highest metal¬ 
licity in our sample, and in turn is the one with the most currently 
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Figure 25. Period versus mass (upper panel), period versus eccentricity 
(middle panel) and stellar mass versus stellar metallicity (bottom panel) 
for all confirmed planets with known values for the stellar type, mass and 
metallicity, and for known values of planetary mass, eccentricity and pe¬ 
riod. There is a total of 317 planets around dwarf stars, 76 around sub-giant 
stars, and 67 around giant stars (including the ones in this work). 


detected planets (two found in this work). This follows the notion 
that the higher the metallicity of a star, the higher the probability 
of finding a planet orbiting them. This result was discussed in R15 
and mimics the trend and dis coveries of planets orbit i ng the most 
metal-rich lower-mass st ars dpischer & Valentill2005l : ISousa et al] 
I 2 OI itlienkins et al.ll2013h . 

One special case is that of HD 47536 b. This planet is orbiting 
a very metal poor giant star ([Fe/H]=-0.65 dex), making HD 47536 
the second most metal-poor star known to host a planet (the most 
metal-poor star was thought to be BD-l-20 2457, with [Fe/H]=- 
1.0 dex, but it was recently shown to be an unstable system; see 
iHorner et alj|201^ . Actually, the probability of finding a planet 
around such a metal poor star is almost 0% (according to Figure 


4 and 7 of R14). The confirmation of a planet orbiting HD 47536 is 
a very interesting case to consider in the study of the formation of 
these systems. 
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RAFT 1: New planetary candidates and updated orbits 


Table Al. RV measurements for HD 11977. 


Table A2. RV measurements for HD 47536. 




FEROS 





FEROS 



JD- 

2450000 

RV 

(m/s) 

Error 

(m/s) 

JD- 

2450000 

RV 

(m/s) 

Error 

(m/s) 

JD- 

2450000 

RV 

(m/s) 

Error 

(m/s) 

JD- 

2450000 

RV 

(m/s) 

Error 

(m/s) 

3100.480 

-54.5 

7.5 

3585.911 

4.9 

8.2 

3454.506 

-65.8 

7.7 

4217.484 

-99.0 

6.9 

3101.488 

-32.4 

8.6 

3586.743 

15.2 

8.7 

3454.509 

-63.9 

8.7 

4231.459 

-125.2 

7.9 

3102.481 

-82.6 

6.9 

3586.903 

-0.4 

7.9 

3460.566 

-39.0 

5.5 

4233.471 

-24.1 

7.0 

3102.486 

-74.4 

7.1 

3590.726 

-5.1 

9.0 

3460.569 

-40.2 

5.3 

4402.754 

150.7 

6.0 

3216.910 

-54.8 

7.5 

3591.767 

0.4 

8.3 

3685.756 

27.3 

8.0 

4407.858 

109.3 

5.6 

3217.842 

-66.7 

11.1 

3591.931 

1.4 

7.2 

3687.752 

52.8 

7.6 

4492.534 

115.9 

10.0 

3218.933 

-59.6 

8.1 

3592.733 

10.0 

9.9 

3688.716 

51.3 

7.1 

4493.533 

110.1 

7.5 

3331.581 

-12.9 

5.9 

3592.933 

3.6 

9.5 

3689.708 

75.0 

7.6 

4494.540 

88.4 

9.9 

3332.556 

-2.8 

5.8 

3593.725 

-11.4 

9.9 

3689.857 

70.3 

7.2 

4509.602 

31.7 

9.1 

3334.622 

11.6 

6.5 

3593.929 

-8.4 

9.0 

3690.663 

-1.7 

7.4 

4510.579 

106.2 

9.5 

3335.688 

6.3 

5.4 

3594.715 

-8.0 

9.0 

3809.589 

-142.1 

6.5 

4572.613 

-75.1 

6.2 

3336.691 

16.3 

5.2 

3594.926 

-2.6 

7.9 

3813.559 

-133.0 

6.8 

4582.500 

-140.6 

6.5 

3338.709 

-12.8 

6.0 

3595.853 

2.5 

14.2 

3869.512 

-108.7 

8.1 

4780.881 

103.1 

4.9 

3448.479 

71.6 

7.4 

3932.943 

-20.7 

7.2 

3873.472 

-114.9 

8.8 

5109.833 

-83.5 

7.5 

3449.481 

107.5 

4.1 

4041.589 

44.3 

6.3 

4040.759 

147.3 

6.7 

5167.827 

-138.1 

5.1 

3449.484 

114.9 

4.8 

4112.543 

115.1 

5.1 

4042.774 

139.5 

6.4 

5201.738 

1.5 

6.4 

3581.730 

11.1 

9.4 

4113.550 

116.7 

3.4 

4112.653 

5.1 

7.6 

5255.591 

74.9 

5.0 

3581.925 

10.6 

8.3 

4115.552 

105.9 

2.4 

4113.636 

4.9 

6.2 

5258.694 

119.2 

5.7 

3582.737 

25.1 

9.7 

4313.859 

-65.9 

16.1 

4114.664 

-12.4 

6.7 

5491.768 

93.1 

9.3 

3582.929 

8.0 

9.0 

4315.791 

-85.3 

8.7 

4115.689 

22.7 

7.0 

5494.891 

92.6 

7.1 

3583.726 

24.2 

9.7 

4316.850 

-78.8 

7.7 

4120.765 

-2.8 

8.0 

5521.825 

-5.6 

7.7 

3583.945 

23.5 

9.0 

4316.927 

-45.1 

7.2 

4121.710 

25.6 

7.5 

5523.764 

-44.3 

7.4 

3584.705 

25.5 

8.9 

4317.857 

-90.2 

14.9 

4122.790 

-28.6 

7.0 

5524.789 

18.9 

7.4 

3585.720 

11.1 

9.3 

4319.871 

-66.3 

8.1 

4123.791 

43.8 

6.9 

5525.881 

-41.4 

7.5 







4124.815 

18.0 

7.1 

5583.808 

-110.4 

7.2 






4157.525 

-23.8 

9.5 

5588.777 

-111.3 

7.2 







4162.558 

65.5 

8.0 

5630.650 

-102.9 

8.6 

JD - 

RV 



RV 














4216.484 

31.3 

6.1 

5632.546 

-117.8 

7.7 

2450000 

(m/s) 

(m/s) 

2450000 

(m/s) 

(m/s) 







3334.620 

49.2 

0.2 

4787.602 

120.8 

0.3 



HARPS 



4045.592 

121.1 

0.6 

4891.528 

-10.6 

0.4 

JD- 

RV 

Error 

JD- 

RV 

Error 

4370.721 

-30.1 

0.3 

5083.811 

-1.5 

0.3 

2450000 

(m/s) 

(m/s) 

2450000 

(m/s) 

(m/s) 

4371.671 

-7.5 

0.4 

5084.748 

-2.9 

0.2 







4622.915 

55.1 

0.4 

5085.754 

-8.5 

0.2 

4045.814 

18.4 

0.3 

4622.442 

-93.4 

0.4 

4684.748 

128.8 

0.3 

5086.730 

-2.1 

0.4 

4085.848 

111.1 

0.5 

4623.462 

-25.5 

0.7 

4785.641 

131.2 

0.6 




4188.594 

-70.4 

1.1 

4684.925 

-18.0 

0.3 







4236.489 

-94.4 

1.3 

4685.925 

7.4 

0.4 



CHIRON 



4370.863 

84.4 

0.3 

4786.762 

49.4 

0.4 







4560.593 

-48.2 

0.3 

4788.863 

39.1 

0.3 

JD - 

RV 



RV 














4561.578 

-54.1 

0.3 

4877.701 

135.4 

0.3 

2450000 

(m/s) 

(m/s) 

2450000 

(m/s) 

(m/s) 

4594.467 

-79.4 

0.3 

4891.645 

66.1 

0.3 

6906.820 

-17.7 

5.4 

6965.553 

-6.1 

4.6 

4596.469 

-101.1 

0.5 

4893.659 

98.7 

0.4 

6919.634 

-16.1 

4.0 

6988.552 

10.6 

5.2 






6921.717 

-22.6 

4.7 

6993.525 

19.6 

5.0 
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6964.740 

-3.5 

4.6 

7013.531 
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5.2 

JD- 
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Error 

JD- 

RV 

Error 







2450000 
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2450000 
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2550.388 
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15.0 

2719.400 

240.6 

10.0 







2622.090 

76.7 

10.2 

2999.379 

10.5 

10.0 

APPENDIX B: 
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2642.576 

178.0 

5.0 

3001.086 

29.7 

10.0 


2669.891 

140.9 

15.0 

3047.180 

113.0 

10.0 

2683.549 

255.8 

10.0 

3052.302 

130.2 

10.0 

2717.693 

192.4 

10.0 

3782.980 

-118.3 

10.0 

CHIRON 

JD- 

RV 

Error 

JD- 

RV 

Error 

2450000 

(m/s) 

(m/s) 

2450000 

(m/s) 

(m/s) 

6953.723 

13.8 

6.8 

6979.718 

-18.5 

6.3 

6959.719 

-105.2 

6.8 

6987.636 

-36.6 

5.8 

6967.768 

48.5 

6.2 

7012.707 

60.7 

6.7 

6977.711 

37.2 

6.9 
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Table A3. RV measurements for HD 110014. 


Table A4. RV measurements for HD 122430. 


FEROS 


FEROS 


JD- 

RV 

Error 

JD- 

RV 

Error 

2450000 

(m/s) 

(m/s) 

2450000 

(m/s) 

(m/s) 

4112.805 

-166.6 

10.5 

4316.480 

91.8 

8.2 

4113.784 

-171.2 

8.0 

4492.899 

120.5 

9.4 

4120.857 

-111.7 

7.4 

4495.852 

106.0 

7.3 

4120.861 

-112.5 

7.9 

5260.801 

53.7 

7.2 

4159.903 

-15.0 

8.3 

5262.783 

42.7 

5.0 

4161.809 

-5.6 

9.5 

5309.667 

40.9 

11.4 

4216.540 

-74.1 

7.2 

5356.570 

156.1 

7.8 

4217.765 

-64.0 

8.9 

5586.890 

-13.0 

5.0 

4217.770 

-64.7 

8.8 

5615.805 

-10.1 

9.2 

4309.474 

79.8 

9.4 

5615.895 

1.0 

8.4 

4309.479 

85.6 

9.3 

5630.733 

-22.2 

7.8 

4312.487 

93.6 

9.4 

5672.706 

-134.2 

7.1 

4313.560 

73.9 

11.2 




HARPS 

JD- 

RV 

Error 

JD- 

RV 

Error 

2450000 

(m/s) 

(m/s) 

2450000 

(m/s) 

(m/s) 


4235.759 

38.1 

0.3 

4891.766 

-12.2 

0.2 

4560.834 

364.5 

0.3 

4931.641 

20.2 

0.2 

4594.768 

289.4 

0.4 

5918.873 

-19.2 

0.2 

4595.565 

284.6 

0.3 

5934.864 

-28.9 

0.2 

4596.735 

287.5 

0.4 

5935.844 

-15.5 

0.2 

4622.599 

189.3 

0.3 

5936.861 

-16.4 

0.2 

4623.663 

204.5 

0.4 

5937.861 

-15.6 

0.2 

4684.471 

272.4 

0.2 

5938.860 

-37.8 

0.2 

4685.514 

247.0 

0.3 







Figure Bl. S-index (upper panel), bisector velocity span (middle panel) 
and FWHM of the CCF (lower panel) versus the measured RV (left panel) 
and versus the residual RV for HD 11977, after fitting the 1-planet solution. 
Only the values for the activity indices laying 2.5cr from the mean value of 
the distribution were considered. The correlation coefficients for the plots 
in the left panels are r= 0.11, -0.29 and 0.05, and for the right panels are 
r = 0.3, -0.19 and 0.25 for the S-index, BVS and FWHM versus the residual 
velocity, respectively. Horizontal dashed lines correspond to the mean of the 
distribution for each plot. 
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Figure B2. S-index (upper panel), bisector velocity span (middle panel) and 
FWHM of the CCF (lower panel) versus the measured RV (left panel) and 
versus the residual RV for HD 70573, after removing the rotation of the star. 
Only the values for the activity indices laying 2.5(T from the mean value of 
the distribution were considered. The correlation coefficients for the plots 
in the left panels are r= -0.19, 0.44 and 0.05, and for the right panels are 
r = -0.28, 0.26 and 0.11 for the S-index, BVS and FWHM versus the 
residual velocity, respectively. Horizontal dashed lines correspond to the 
mean of the distribution for each plot. 
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Table A5. RV measurements for HD 70573. 


FEROS 


JD- 

2450000 

RV 

(m/s) 

Error 

(m/s) 

JD- 

2450000 

RV 

(m/s) 

Error 

(m/s) 

3099.524 

-45.7 

11.7 

4114.743 

17.1 

11.0 

3101.539 

-70.7 

10.3 

4114.835 

39.8 

9.9 

3449.674 

29.2 

16.4 

4115.721 

-30.3 

49.1 

3687.838 

5.3 

11.2 

4121.748 

90.1 

8.8 

3807.576 

-48.4 

12.8 

4121.856 

106.3 

12.9 

3807.638 

-77.8 

23.7 

4122.802 

82.6 

11.7 

3809.735 

-9.4 

22.8 

4122.849 

77.9 

8.5 

3810.604 

-61.5 

12.3 

4123.816 

-19.8 

10.6 

3813.660 

-39.7 

11.3 

4123.859 

-11.1 

13.4 

3816.663 

-16.1 

11.8 

4124.790 

19.0 

10.1 

3817.716 

-75.1 

16.1 

4124.866 

26.8 

10.6 

3869.524 

14.4 

10.4 

4159.579 

-10.1 

16.3 

3871.514 

4.4 

12.7 

4159.769 

13.5 

21.7 

3872.493 

-23.1 

9.9 

4160.718 

-26.2 

15.7 

3873.511 

-113.9 

16.1 

4162.588 

-14.8 

9.8 

4039.874 

79.8 

10.1 
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-16.1 

14.2 

4040.875 

-51.3 

7.2 

4163.601 

-71.2 

17.1 

4041.799 

42.6 

12.5 

4164.550 

51.0 

11.5 

4042.873 

-20.1 

5.1 

4164.691 

53.0 

15.4 

4043.863 

51.1 

9.4 

4165.552 

3.1 

11.4 

4048.838 

48.3 

8.4 

4165.683 

-24.0 

15.2 

4053.864 

-11.7 

10.2 

4216.511 

-52.2 

8.6 

4112.688 

67.7 

8.7 

4216.571 

-61.6 

12.3 

4112.783 

72.0 

12.2 

4217.527 

26.7 

10.3 

4113.652 

-11.0 

11.8 

4407.787 

14.8 

9.8 

4113.734 

-31.6 

12.9 

4412.843 

69.3 

12.0 

4113.839 

-83.8 

11.6 

4413.747 

-20.2 

16.2 

4114.645 

4.2 

13.2 

4504.745 

38.6 

17.3 



Figure B3. S-index (upper panel), bisector velocity span (middle panel) and 
FWHM of the CCF (lower panel) versus the measured RV (left panel) and 
versus the residual RV for HD 122430, after fitting the best 1-planet solution 
we were able to obtain. Only the values for the activity indices laying 2.5(T 
from the mean value of the distribution were considered. The correlation 
coefficients for the plots in the left panels are r= -0.17, 0.01 and 0.04, and 
for the right panels are r = -0.32, 0.1 and 0.04 for the S-index, BVS and 
FWHM versus the residual velocity, respectively. Horizontal dashed lines 
correspond to the mean of the distribution for each plot. 
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